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Oh. Abstract: We study the value of shadowing corrections (SC) in HERA kinematic region 

r-| ! in Glauber - Mueller approach. Since the Glauber - Mueller approach was proven in pertur- 

bative QCD in the double logarithmic approximation (DLA), we develop the DLA approach 
for deep inelastic structure function which takes into account the SC. Our estimates show 
small SC for F2 in HERA kinematic region while they turn out to be sizable for the gluon 
structure function. We compare our estimates with those for gluon distribution in leading 
order (LO) and next to leading order (NLO) in the DGLAP evolution equations. 
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In this letter, we investigate the role of the shadowing corrections (SC) on the value 
of the deep inelastic structure functions and on the scaling violation mechanism in HERA 
kinematic region. Our estimates of the SC are done using the Glauber - Mueller formula. 
This formula was proven by Mueller [|l| in the DLA of pQCD and was studied in details in 
Refs.0|]^[^]. The structure function F2{x,Q'^) reads as: 



where x is the Bjorken scaling variable and is the photon virtuality. is the colour 
number and Zj is the charge fraction of each quark, Nj is the number of flavours taken into 
account in the quark- nucleon scattering and 6^ is the impact parameter for the scattering 
of the quark - antiquark pair with transverse splitting r±_ with the nucleon target. 

The opacity function VL is, generally speaking, an arbitrary real function and it has a 
simple physical interpretation: is the probability that the quark-antiquark pair do not 
suffer an inelastic scattering. This function should be determined in QCD. For the DLA it 
was shown ( see Refs.[Q [Q) that the impact parameter dependence can be factorized out 
and Vt can be written 

A 2 

^{x, Q\ h^) = -^xG{x, Q')S{b^) (2) 

where S{b±) is the profile function for a nucleon with radius R, and is taken as an exponential 
function, S{b±) = -:^exp{—h\/I^), in our calculations. 

It should be stressed that Eq. (|^) takes into account the unitarity constraint and with 
VL defined in Eq. (Q) gives the DGLAP evolution equation in the kinematic region where 
< 1. 

In order to investigate the x and evolution of the gluon distribution in this approach, 
we will consider first our Born term with the gluon distribution taken in the DLA limit. 
Since the Born term is equivalent to the DGLAP expression for F2 in the DLA limit, we will 
take xG^^^ as in Ref.|^ and Thus, the gluon distribution reads 

xG^^^(x,g^) = Go/o(l/), (3) 

where the variables are y = 27^/n(l/x)/n(t/to) and t = Ini^Q"^ / h?). The QCD constants are 

7 = y'l2//5o, /5o = 11 - ^ and A = 0.232 GeV^ . The constant Gq plays the role of the fiat 
initial condition, since the Bessel function Io{y) goes to 1 as ?/ goes to zero. We disregard 



2 



the sub-leading corrections to the DLA gluon proportional to a^lnQ"^ proposed in Ref.[0. 
Integrating expression (|l]) over the Born term then reads 

F. = F,{x, Ql) + A ^as(g'^)xG^^^(x, Q'') , (4) 

where we have taken Nf = 3. The expression (^) gives the sea component of F2 generated 
by the gluon evolution from the initial virtuality Qq to Q^. Going from expression (|l|) to 
we have taken Q'^ = Thus, the lower limit works as a cut off for the large distances 

effects over F2. These effects are included in F2{x,Ql), the value of the structure function 
for the virtuality Q^. It has a nonperturbative origin and takes into account the amount 
of qq pairs not generated by the perturbative transition g —>■ qq. We will parameterize the 
initial structure function by the expression 

F2{x,Ql)=CoX-'''{l-xY\ (5) 

where Co is a constant that adjusts the nonperturbative contribution. This expression re- 
produces the soft pomeron behaviour as x — > 00) presented by the 7* — Nucleon 
cross section in the low Q"^ region 0. Since we have used A^/ = 3, we should add the charm 
component . This component is generated perturbatively from the 7* — gluon fusion 
mechanism with the gluon distribution given by the DLA expression (|^). This mechanism 
is discussed in detail in Ref . [|T0| . Finally, we obtain the following expression for F2 

F2{x, Q') = F2{x, Ql) + F^;^l^ + FIj,^^ . (6) 

To fit the expression (Q) to the HERA data, we have taken the F2 points which lies in the 
region IGeV^ < < 100 GeV"^ and x < 10~^, where we expect that our DLA approach 
to SC is valid. The HI and ZEUS results were taken from Refs.[jll| and [|12|, respectively. 



In figure (|I|) we present the fit for a subset of the experimental data. The parameters 
used are Go = 0.136, Cq = 0.273 and = 0.330 GeV^^. The values of the parameters were 
chosen in such a way to minimize the x^, which corresponds to /d.o.f. = 124/222. We can 
see from the figure that the steep behaviour of the deep inelastic struture function is well 
described by the DLA evolution of the gluon distribution, regarded we have included enough 
nonperturbative qq pairs. With this set of parameters also the scaling violation of F2 
can be described, as shown in figure (0). Taking a small value for the initial virtuality we 
can generate the DLA behavior for ^ 2GeV^. A similar result was obtained in Ref.p], 
but in a not completely DLA limit]^. It is important to note that our aim in this letter is to 
describe HERA data in a completely consistent DLA limit, and not to provide an overall fit 
to existing high energy data. 

*The authors have taken Pqg — + {1 — z)^ and the sub-leading factor (t/to)^^ with S = (11 + ^^)//3o 
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Since we have described the data with the Bom term of the DLA expression (|I]), we 
can investigate the amount of shadowing corrections predicted for F2 on Glauber - Mueller 
formula of Eq. (|I]). For that, we substitute the Born term in expression (^) by the full 
series, which is taken into account in expression (|lD. In figure (J^) we present the results for 
F2 as a function x, and in figure (j^), as a function of Q^. As we can see, the shadowing 
corrections are important only for very small values of x and moderate values of Q^. We 
would like to recall that in Eq. (|T]) we put the upper limit of integration equal to or in 
other words we consider only the SC which are originated from sufficiently short distance, 
namely, r± < 1/Qo ~ 0.35/m. In fact, we do not take into account the SC at large distances 
considering that they have been included in the initial parton distribution of Eq. (^. 

Therefore, we are calculating only perturbative shadowing. In the kinematic region of 
present data, the corrections lie within the experimental error. 

We plot also in figures and (|^) the SC for F2 predicted by the Glauber approach taking 
into account the leading order (LO) and next to the leading order (NLO) gluon. In both 
cases, we have used the modified Mueller formula discussed in Ref.0. In this formula, the 
Born term is taken in leading aslnQ^ aproximation (LLA(Q^)), while the correction term is 
taken in DLA. For practical purpose, we use the structure function F2, solution of DGLAP 
equations, as the Born term in expression (^). We have taken only the GRV distribution 
since those distributions evolve from small virtualities and can be compared with our DLA 
approach. We see from the figures that the LO gluon predicts much more SC to F2. It means 
that the scaling violation suffers a stronger modification for the LO gluon when compared 
to the simple DLA gluon and to the NLO gluon. 

The Glauber - Mueller approach cannot be considered as a full description of the SC, 
because it was assumed that only quark - antiquark pair embodies a multi rescatterings with 
the target. As was shown in Refs. 010 gluon rescatterings turn out to be more 
essential. To demonstrate this fact we calculate here the Born term of Eq. ([l|) but given by 
expression: 

F2 = F2{x, Ql) + A ^as(Q'^)xG^^(x, Q'') (7) 

where xG"^*^ is the gluon structure function calculated in the Glauber - Mueller approach, 
namely 

xG^"^ x,Q2 = % ^ dhl{l - e ^ } 8 

TT^ J-L r I Jx x' Jo 

where the opacity VLq = Expression (H) is the Mueller formula which was discussed 

in detail in Ref. [Q. When Eq. is included in expression (|^, the Born term reproduces 
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Eq. (^), since the Born term is the DGLAP equation in the DLA hmit. The other terms take 
into account the shadowing corrections to the gluon distribution. The results are shown in 
figure (|^). Comparing figures and one can see that the SC due to gluon rescattering 
is bigger then the corrections due to quark rescattering. 

In order to complete our discussion, we plot in figure (H) the DLA gluon distribution 
given by expression and the corrected gluon distribution given by the modified Mueller 
formula The LO and NLO gluon distribution given by the parameterization GRV95 are 
ploted also. As we can see, the DLA distribution predicts an amount of gluons closer to the 
NLO DGLAP evolution. It is not a coincidence, since the NLO GRV distribution has a fiat 
behaviour for = OAGeV^ while the LO gluon distribution has already a steep behaviour 
in the small x region for this low value of Q^. 

Fig. (P) shows the main conclusion of this letter: the SC turns out to be big ( about 
40% - 50% ) in the gluon structure function but their manifestation in F2 is rather small as 
we have discussed ( see Fig. (|]) ). Comparing also figures (0), (|[) and we can see that 
the SC for F2 have a strong dependence on the ammount of gluons taken into account in the 
QCD evolution. This conclusion calls for new measurements in the high energy kinematic 
region more sensitive to the value of the gluon structure function than the measurements of 
F2. 

ALAF- acknowledges CAPES and MBGD acknowledges CNPq for partial financing. 
EML thanks E. Gotsman and U. Maor for everyday discussions on the subject. 
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Figure 1: The structure function F2 evolution from the scahng violation mechanism in DLA 
aproximation as a function of x. The solid line represents the Born term and the dashed line 
includes the shadowing corrections (SC). 
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Figure 2: F2 evolution from the scaling violation mechanism in DLA as a function of 
( for the scaling violation figures, the value of i goes from for x = 2.0 x 10~^ to 7 for 
a; = 8.0 X lO'^). 
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Figure 3: The structure function F2 evolution in LO and NLO as a function of x. The Born 
term (DGLAP evolution) for F2 NLO and LO numerically coincide for Q"^ > lOGeV"^. 
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Figure 4: F2 evolution from the scaling violation mechanism in LO and NLO as a function 
of Q\ 
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Figure 5: F2 evolution from the scaling violation mechanism in DLA with the gluon distri- 
bution xG^^ given by the Glauber - Mueller approach. 




Figure 6: The gluon distributions in DLA with and without SC - calculated from the 
Glauber - Mueller approach, compared with LO and NLO gluon distribution from GRV95 
set. 
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